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Abstract We analyzed an XMM-Newton EPIC observation of a bright knot, named FilD, in the northern rim of the Vela SNR, 
where the shock has encountered a cloud. The good combination of sensitivity, spectral, and spatial resolution allowed us to 
describe the internal structure of the observed ISM clouds and to obtain estimates of their temperature, density, O, Ne, and Fe 
abundances, and of their extension along the line of sight. We also examined the interaction of the shock with the FilD knot and 
estimated that the time elapsed from the shock impact is about one cloud crushing time. Our analysis allowed us to conclude 
that the observed X-ray emission is best explained by the propagation of transmitted shocks through ISM inhomogeneities with 
an inward increasing density profile. We found that the interaction with the shock determines the heating of the inner part of 
the clouds and the evaporation of the outer part. 

Key words. X-rays: ISM - ISM: supernova remnants - ISM: clouds - ISM: kinematics and dynamics - ISM: individual object: 
Vela SNR 



1. Introduction 

Supernova remnants (SNRs) play a leading role in settling the 

■ mass distribution and the energ y exchanges in the i nterstellar 
medium (ISM) and, according to lCox & Smithl ( fl974 ). they ac- 
count for the heating of about 1/10 of the volume of the inter- 
stellar medium of our galaxy. During their expansion, SNRs 
heat, compress, and ionize the ISM and determine the destruc- 

^ tion of the interstellar dust grains and the consequent release of 

■ heavy elements in the interstellar gas. 

The inhomogeneities of the ambient medium influence the 
direction and the speed of the SNRs expansion. The optical and 
X-ray emissions of middle-aged SNRs arise from the interac- 
tion between the shock and the ISM clouds: optical filaments 
arise from slow radiative shocks propagating in a dense (n ~ 10 
cm""*) and cool (T ~ 10^ - 10^ K) medium, while X-ray emis- 
sion originates from hotter and less dense regions. However, the 
details of the heating processes are still not well understood. 

Three diff'erent physical scenarios have been proposed to 
explain the clumpy emission of the SNRs. A first scenario 
associates the optical emission to a transmitted shock travel- 
ling through dense clouds and explains the X-ray emission 
with the evaporation of the inhomogeneities engulfed by the 
main blast wave. This model has been applied in the past to 
observations of SNRs ( Fesen etaL..1982. Chai'les et aL..1985l 
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ICraig et alJll997h an d, more recently, to describe the compos- 
ite SNR G327.1-1.1 jBocchino & Bandieral2003l) . 

A second class of models points out the importance of 
the overpressured regions in the post-shock flows and asso- 
ciates the X-ray emission to the inter-cloud medium heated 
by reflected shocks. Hydrodynamical simulations show that 
the cloud density contrast drives reverse reflected shocks 
in the inter-cloud medium (see e.g. Stone & Normaril ll992t) . 
According to this model the X-ray emission is localized "be- 
hind" the optical emission, which is associated with the shock 
transmitted into the inhomogeneities. This scenario has been 
extensively adopted to explain the optical and X-rav enussion 
of the Cvgnus Loop ('Gra ham et alJl995llLevenson et alJl996l 
iMivat a & Tsunemi 2001). 

A third interpretation is based on the propagation of trans- 
mitted shocks through inhomogeneities with an inwards in- 
creasing density profile. These shocks can account for the ob- 
served X-ray emission when the cloud contrast is low and/or 
when the prim arv blast wave shock is fast. This scenario was 
proposed bv .Bocchino et al.l J2OOO) to explain the optical and 
X-ray emission of two regions of the Vela SNR in which the 
main blast wave seems to interact with a small isolated cloud 
(FilD region) and with a cavitv wall (FilE r egio n), and has been 
recently applied by Patnaude et alJ ( l2002i) and lLevenson et all 
( l2003t) for two regions of the Cygnus Loop. 

The instrumental performances of past generation X-ray 
satelhtes did not always allow discrimination between these 
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models or an unambiguous description of the emission mecha- 
nism in middle-aged SNRs. 

The Vela SNR is the nearest object of its kind, so it rep- 
resents a privileged laboratory for the study of the shock- 
cloud interaction. The distance of the Vela SNR is ~ 280 pc 
J " < I ■ 

JBocchino et al. 1999, Cha et al. 1999) and at this distance one 
arcmin corresponds to only 0.08 pc. The age of this SNR has 
been estimated to be ~ 10000 yr, in good agreeme nt with 
the characteristic age (~ 11200 vr. iTavlor et al.lfl993l) of the 
pulsar PSR B083 3-45, w hich is associated with the remnant 
( Weiler & Panagia ll980l and Weiler & Sramek 1988). 

The northern rim of the Vela SNR presents all the character- 
istic features of a complex interaction between the blast wave 
shock and several inhomogeneities of the ambient medium, as 

S iointed out b y the ROSAT observations of FilD and of FilE 
Bocchino et al 1994. .1997.. ,1999, 2000) . Analysis of the X- 

ray emission of these regions showed that the observed spec- 
tra are described significantly better by a 2-component opti- 
cally thin thermal plasma in equilibrium of ionization than by 
a single-temperature model or a non-equilibrium of ionization 
(NEI) model. As for FilD, the two components with temper- 
atures 1.10 + 0.06 X 10^ K and 6+3 x 10'' have been associ- 
ated with the shocked plasma of the cloud and with the inter- 
cloud medium, respectively. However the instrumental charac- 
teristics of the PSPC detector did not allow us to describe the 
shock-cloud interaction univocally or to get detailed informa- 
tion about the morphology of the emitting plasma and about its 
chemical abundances. 

Here we present analysis of an XMM-Newton EPIC obser- 
vation of the FilD region. The aim of this work is to exploit the 
EPIC characteristics to derive a better description of the shock- 
cloud interaction process in this region of the Vela SNR and to 
discriminate among the different models described above. 

The paper is organized as follows: in Sect.|2]we present the 
EPIC data analyzed here, while in Sect.|3lwe describe the mor- 
phology of the X-ray emission, identifying the physically dif- 
ferent emitting regions, and then present the results of the spa- 
tially resolved spectral analysis. In Sect.l^we study the thermal 
structure of the source, show the 3D map of the plasma, and 
discuss the shock-cloud interaction in this region of the Vela 
SNR. The presence of NEI effects is examined in Sect. |4] and, 
finally, our conclusions are drawn in Sect.ls] 

2. The Data 

The XMM-Newton observation of Vela FilD we present is a 
Guaranteed Time O bservation, performed with t he EPIC MPS 
/Turner et al J 12001 ^ cameras and the EPIC pn (Striider et al.' 
I2OOI ) camera on 9 November 2001 and with pointing coor- 
dinates a (2000) = 8''35'"44', 5 (2000) = -42°35'29". EPIC 
cameras are CCD detectors with spectral resolution E/AE ~ 
20 - 50 in the 0. 1 - 10 keV band. The field of view (FOV) has 
a diameter of ~ 30', and the pixel size is 1.1" (MOS) and 4.1" 
(pn), while the mirror PSF width is 6" - 15" FWHM-HEW. 
The XMM-Newton mirror assembly is subject to vignetting, so 
there is a decrease in effective area with increasing in off-axis 
angle. 

Table [1] summarizes the main information about the data. 



Table 1. Relevant information about the data. 



CAMERA 


f„p (ksec)' 


Mode 


Filter 


Ki 


pn 


26.8/11.5 


Ext. Full Frame 


medium 


291115 


MOSl 


31.1/25.2 


Large Window 


medium 


203113 


M0S2 


31.1/25.2 


Large Window 


medium 


207807 



I. Unscreened/Screened exposure time. 



II. Number of events in the whole field of view (0.3 - 10 keV). 



2.1. Data processing 

Data have been processed with the Standard Analysis System 
(SAS V5.4. 1) pipeline. We further screened the data by exclud- 
ing some residual hot pixels (especially from the MOSl ob- 
servation) and by removing the time intervals which showed 
anomalous background enhancements (probably due to soft 
protons). In particular, we extracted the background light curve 
with a bin size of 10 s, at energy > 10 keV and we eliminated 
time intervals with high counts (> 30 for pn data and > 15 
for MOSl and MOS2 data). Table[2shows the unscreened and 
screened exposure time tg ^p^, the pn observation is strongly con- 
taminated by proton flares (~ 57 % of the total exposure). 

Spectral analysis was performed using XSPEC. The an- 
cillary response files (arf) were generated with the SAS 
ARFGEN task, and t he event files were processed using the 
EVIGW EIGHT task dArnaud etalJ200l l) to coiTect vignetting 
effects. We summed the MOSl and MOS2 spectra averag- 
ing the arf files (and summing the exposure times) and using 
the MOS 1 response file for the analysis of the summed spec- 
trum. The fittings were performed simultaneously on the pn and 
MOS-sum med spectra. All the r eported errors, calculated ac- 
cording to lLampton et a D lll976l) . are at 90% confidence level. 

2.2. Background subtraction 

We removed the X-ray sources present in the EPIC field 
of view and not correlated with the Vela SNR. Detection 
of these sources was performed on the summed data set 
(MOSl-HMOS2-i -pn), using the Wavelet Transform algorithm 
PWXDETECT llDamiani et alJl997albl) . In this way we found 
position, extension, and X-ray flux of 21 sources; 10 among 
them have an optical counterpart in the USNO-A2.0 archive. 
All the detected sources have a significant signal only in the 
energy range 1-10 keV, so we subtracted their contribution 
in this band. Since the Vela emission covers the entire FOV, it 
was not possible to isolate a region where the local background 
could be extracted. 

The contribution of the incident photon X-ray background 
was derived from the EPIC public observation 0106260201, 
whose FOV is wholly outside the Vela SNR shell, and pointed 
toward the source RXJ0806.4-4123, which is only ~ 5° distant 
from our observation. In the 0.3 - 2 keV band (where our anal- 
ysis is focussed) the count rate contribution of this background 
is between ~ 5 % and ~ 15 % of the total signal. In particular. 
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in a 4" bin, the mean 0.3-2 keV MOS count rate of the photon 
background is R ~ 8.6 x 10"'' s"' with a mean photon energy 
E ~ 670 eV; we used these values to compute the background 
subtracted E vs. R plot described in Sect. 13. 21 We also used this 
observation to extract a pn and a MOS background spectrum 
for the spectral analysis. However, as shown in Appendix A, 
the best fit parameters are not significantly influenced by the 
choice of the background. 

As for the non-photon background, which is probably asso- 
ciated with a residual flux of soft protons, we found that above 
1 keV the non-photon background amounts to > 95% of the 
total signal for the pn camera and > 70%, for the two MOS 
cameras and that: i) In the 0.3 - 0.5 keV band the non-photon 
background contribution is too low to be estimated; ii) in the 
0.5 - 1 keV band it amounts to 2 - 12 % of the total signal for 
the pn camera, and iii) in the 0.3 - 2 keV band it is negligible. 



3. The data analysis 

3.1. The region morphology 

To reduce statistical fluctuations in the images, we imple- 
mented a weighted average of the pn, MOSl, and MOS2 count 
rate maps in three different bands: 0.3 - 0.5 keV, 0.5 - 1 keV, 
and 0.3 - 2 keV. The first two bands were chosen to understand 
how the emission is distributed in the FOV at different energies, 
while the image in the 0.3 - 2 keV band allowed us to make 
a comparison between the EPIC and the PSPC observations. 
All images are MOS -equivalent, and the weight of the different 
cameras were chosen taking the thermal spectrum of the source 
into account. We estimated that in all these images the relative 
statistical error of the count rate in each bin is < 18%. 

Fig. [0 shows the count rate images (vignetting corrected) 
in the 0.3 - 0.5 keV band and 0.5 - 1 keV band; these im- 
ages are adaptively smoothed with cr„„„ = 4" and cr,,,^. - 24". 
North is up and E ast is on the left. Not ice that, assuming a 
distance of 280 pc jBocchino et aLlfTQQOl) . the diameter of the 
EPIC FOV corresponds to an extension of ~ 2.4 pc. There are 
striking differences in the morphology of the emission in these 
two bands. The FilD complex, in the middle of the FOV, is 
clearly visible in the soft band {left panel Fig. [0, and its in- 
ternal structure is resolved in good detail with XMM-Newton. 
FilD looks composed by two bright parts (A and B in Fig.[0. 
South of FilD there is a wide region, labelled C, that extends 
westwards. As shown in the ROSAT observation RP500015 
llBocchino et al.lll999t Fig. 2b), region C stretches beyond the 
EPIC FOV and is probably part of a very large structure. FilD 
is instead barely visible in the hard band {right panel), where 
another large structure is present. This bright structure, named 
RegNE, has a slightly higher count rate than FilD in the 0.3 - 2 
keV band (Fig.|2} and was not visible in the ROSAT PSPC ob- 
servation. 

It is interesting to note how EPIC data relate to the op- 
tical observations of the FilD region. Figure |3] shows a color 
composite optical image of FilD (from lBocchino et"al Ifcooo), 
where we have superimposed the X-ray contour levels of the 
EPIC observation in the 0.3 - 0.5 keV band. The FOV of this 



observation is indicated by the white square in the left panel of 
Fig.E 

In the southwestern part of the image there is a clear rela- 
tionship between optical and X-ray emission: maximum X-ray 
emission contour levels, in fact, are just outside the optical fila- 
ment which follows a path that "grazes", without crossing, the 
two FilD regions with the highest X-ray surface brightness (A 
and B in Fig.Q. This kind of rela tionship between optic al and 
X-ray emission, also observed bv lPatnaude et alJ (l2002h in the 
southwestern region of the Cygnus Loop, suggests that a trans- 
mitted shock is propagating through a cloud with an inwardly 
increasing density profile, generating optical emission in the 
inner and denser part of the cloud and X-ray emission in the 
outer part. 

3.2. Analysis of the photon energy map 

The XMM-Newton data allowed us to study the distribution of 
photon energies across the FOV and to produce a high resolu- 
tion mean photon energy {E) map, shown in Fig.|2l right panel. 
To realize this map, we considered the MOS 1 and MOS2 event 
files in the 0.3-2 keV band; a bin size of 10" allowed us to col- 
lect more than 10 counts per pixel everywhere. For each pixel, 
we calculated the mean photon energy and then we smoothed 
the map with a gaussian distribution with cr = 20" (the error 
on each value of £ is ~ 10 eV). The map clearly shows that E 
has a minimum in the FilD region and a maximum in RegNE. 

The large inhomogeneities in the count rate and mean en- 
ergy maps suggest that the blast wave shock propagates in a 
complex environment, generating different thermal conditions 
in the shocked plasma. To get information about the distribu- 
tion of the physical parameters of the plasma and about the ge- 
ometrical morphology of the inhomogeneities, we investigated 
the relation between the mean photon energy in a given pixel 
and the corresponding count rate R. 

Figure |4] shows the observed E vs. R relation (background 
subtracted), where each point corresponds to a 10" pixel of the 
mean energy map. The global trend is a decreasin g one, as 
expec ted in the case of pressure equilibrium ( B occhino et alJ 
l2004 . It is interesting to note how the different structures 
present in the FOV are disposed in the graph. The black points 
correspond to the photons emitted from FilD and an adjacent 
region; they clearly identify a very narrow strip, suggesting that 
FilD and the surrounding fainter regions are part of a single 
structure. We can then assume that FilD is a largely isobaric 
structure formed by components with different temperatures 
and/or extensions along the line of sight. The red points, cor- 
responding to the southeastern region, distinctly identify a dif- 
ferent flatter trend, with a higher E, hence a higher mean tem- 
perature along the line of sight and a higher value of pressure 
and/or volume cBocchino et al. 2004.). This region may there- 
fore be physically different from FilD. Finally, the events de- 
tected in the northeastern region of the FOV (blue) have the 
highest mean energy and the decreasing trend is less marked. 
These points do not overlap with the points either of FilD or 
of the southern region, suggesting that we are observing three 
physically different ISM inhomogeneities. 
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Note that a plasma heated by a reflected shock would yield 
an increasing trend in the E vs. R scatter-plot. In fact, the hot- 
ter plasma behind the reflected shock front is also more com- 
pressed since it was shocked twice. Therefore, in the reflected 
shock case, the plasma with the higher temperature should also 
have a higher emission measure. Thus the E vs. R trend, visible 
in Fig. 13 suggests that we are observing ISM clouds heated by 
a transmitted shock. 

3.3. Spectral analysis 

The analysis above allows us to define spatial regions in the 
FOV with homogeneous physical properties across the plane 
of the image. We can now perform a spatially resolved spectral 
analysis on them to obtain the physical parameters. We ana- 
lyzed the spectra extracted from the 16 subregions indicated in 
the right panel of Fig. 121 These subregions cover all the signif- 
icant structures in the FOV, identified on the basis of the mean 
photon energy map of Sect. 13.21 and their shape and size were 
chosen so as to have enough photon statistics (> 2000 counts 
per subregion) and very small fluctuations of the mean photon 
energy (< 4%), which is almost uniform in each subregion. In 
this way we can probe the different spectral features between 
regions whose extension is only a few arcminutes and which 
do not have internal spectral inhomogeneities across the plane 
of the image. Subregions 1-6 are located on FilD; subregions 
7-12 cover the south western structure; subregion 13 is in 
the "transition zone" between FilD and RegNE which is cov- 
ered by subregions 14 and 15; subregion 16 corresponds to the 
northern area with very low surface brightness. 

All the extracted spectra are described well by two 
MEKAL components of an optically-thi n thermal plasma 
in coUisional ionization equilibrium, CIE dMewe et al.lll985l 
iMewe et al.lll986l iLiedahl et alJIlQQSb . with non-solar Ne and 
Fe a bundances and modified by the interstellar absorption of 
iMo rrison & McCammon (1983). Figure [S] shows a represen- 
tative spectrum, extracted from subregion 7, with its best fit 
model. 

The elements whose Unes mainly contribute to the observed 
spectra are O, Ne, and Fe. In particular in all spectra the K-shell 
Une complexes of O VII (at 0.56 keV), O VIII (at 0.65 keV), 
and Ne IX (at 0.92 keV) are visible. As for the O abundance, 
we found in each subregion 0/Oq = 0.9 -H 1.1 so we tied this 
abundance to 1 . The blending between the Fe XVII L-lines and 
the Ne IX lines generates an entanglement between the best fit 
values of the abundances of Ne and Fe, and the statistics were 
not enough to discriminate between these two values in most 
of the regions. To overcome this problem we focussed our at- 
tention on the subregions 14, 15, and 16 where the signal to 
noise ratio is best at energies of ~ 0.9 keV. In these subregions 
the mean value of the ratio (Ne /NeQ)/{Fe /FeQ) is 4.4 and we 
assumed that this ratio remains the same in the other subre- 
gions. Therefore, in all the subregions, we left the model Fe 
abundance free and linked the Ne abundance to it, so as to have 
{NelNe^)l{FelFe^) = 4.4. 

It was not possible to obtain a good estimate of the Nh pa- 
rameter (we found Nh - l-Sl] 5 x 10^" cm"^), so it was fixed 



to 1 X 10^" cm ^, which is a typical value in this part of the 
Vela SNR (Bocchino et al. 1994, Lu & Aschenbach 2000, and 
iKahn etalJl985l). 

Figures |6}-|8]describe in detail the findings of our analysis. 
The achieved results can be summarized as follows: 

1. In the subregions 14 and 15, spectra can be modelled with 
one component. 

2. The temperatures Tj and Tu of the two components are 
both rather uniform in the FOV (Fig.0 with mean values 
(between the subregions): T/ = 1.14 + 0.02 x lO*" K (cooler 
component) and Tu - 3. 0+0. lxlO^K(hotter component). 

3. The emission measures of both components have large vari- 
ations, as shown in Fig.|8] 

4. A third MEKAL component in the model did not improve 
the description of the data. 

3.4. We and Fe abundances 

Figure |6] shows the best-fit values of Ne/Ne^ in the 16 subre- 
gions versus the mean MOS count rate (for a bin size of 10") 
in the corresponding subregion. The Fe abundances can be de- 
duced from the condition (Ne/NeQ)/{Fe/FeQ) = 4.4. 

In 10 subregions we found Ne abundances significantly 
higher than cosmic value, with an average Ne/NeQ = 1 .7 + 0.2 
over all the subregions. Typically the Ne overabundance is as- 
sociated with overabundance s of O, Mg, Si (Sasaki et al. 200ll) 
and Fe jMivata et alJ I2OOII) and is a clear signature of the 
ejecta, but this is not the case, since Fe is underabundant and 
we found solar O abundances. Similar values of the Ne abun- 
dance in the Vela SNR were found bv lPlucinskv et al] l2002l) . 
who showed that the Vela Bullet D X-ray emission can be 
described with a non-equilibrium of ionization model (which 
seems to exclude that Bullet D is a fragment of ejecta) with 
Ne/Nen » 2. Moreover, in a small region of the Cygnus Loop, 
Mivata & Tsunemi 2001 found analogous O, Ne, and Fe abun- 
dances (Ne/Nee = 1.7!|J-^, 0/Oe = 0.8!°'^ and Fe/Fee = 

In all the analyzed subregions the Fe is underabundant 
with an average va lue Fe/Fe^ = 0.39 + 0.05. According to 
IVancura et al. ( 1994) (VRD), the fraction of Fe locked up in the 
pre-shock ISM grains is 0.98, so in the X-ray emitting plasma 
at the shock front we have Fe/Fee = 0.02. VRD used a nu- 
merical model to describe the evolution of the grain destruc- 
tion which generates the release of Fe in the interstellar gas in 
the post-shock region. However, assuming the typical values of 
this zone of the Vela shell, shock inter-cloud velocity v,c = 700 
km/s, and inter-cloud density n,v = 0.05 cm""* (Sect. l4.4"Tl . this 
model predicts that the observed Fe abundance can be obtained 
only after ~ 20000 yr, i. e. about two times the age of the Vela 
SNR. 

In conclusion, it is difficult to explain the observed Ne and 
Fe abundances, and it is necessary to obtain more information 
about the value and the spatial distribution of the abundances 
of these elements in the Vela SNR. 
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4. The physical interpretation 

4.1. Thermal structure of the region 

As we pointed out in the previous section, the observed spec- 
tra are described well by two thermal components at ~ 10^ K 
and ~ 3 X 10^ K. Since in each subregion the mean photon 
energy is almost uniform (maximum variations of ~ 25 - 30 
eV), two distinct components necessarily indicates two differ- 
ent thermal phases of the interstellar medium intercepted along 
the line of sight. Since the emission measures per unit area are 
highly inhomogeneous, we can conclude that the EPIC data 
reveal for the first time that both these phases are associated 
with the ISM inhomogeneities (therefore, in what follows, we 
will use the term "clouds" to indicate these hot. X-ray emit- 
ting struct ures). This is at varia nce with ROSAT data, which 

— — II 

suggested ( Bocchino et al. 1999) an association of the hotter 

thermal component with the inter-cloud medium (the tenuous 
and almost uniform ISM phase which surrounds the clouds). 
Clouds can so be described as two-temperature structures. 

It is also reasonable to assume that the temperature is 
uniquely d etermined by the local d ensity, because the shock 
equations jLandau & Lifshitj|l959l) predict that the temper- 
ature increase of the shocked plasma is proportional to the 
square of the Mach number M of the shock, which changes 
with the local density n (M oc «-'^^). We conclude that, since 
the temperature of each component is uniform in the FOV, the 
density of each phase is also similar for all the observed clouds. 
Therefore, differences in surface brightness imply differences 
in the extension of the two phases along the line of sight. The 
two components of the X-ray emission can be related to two 
different parts of the ISM clouds with different densities: a core, 
which corresponds to the cooler and denser component, and a 
hotter and less dense corona, which reasonably surrounds the 
core. 

Figures0and|Slshow the best-fit temperatures and emission 
measures versus the mean MOS count rate (for a bin size of 
10") in the corresponding subregion. These results allow us to 
conclude that the trend in the E vs. R scatter plot of Fig.l^is 
essentially determined by the (relative) dimensions of the cores 
and of the coronae. 



4.2. Density of the cores and of the coronae 

An upper limit to the density of the cloud cores (and hence 
of the coronae) is, obviously, the density nopi of the optical 
filament of Fig. |3] T he multi-wavelength s tudy of the FilD 
region performed by iBocchino et alJ (l2000l) has shown that 
nop, - 3.2-8.2 cm"^, so we infer that the density of the X- 
ray emitting plasma must be $ 5 cm"^. 

Considering the lower limits (at 90% confidence level) of 
the emission measure (shown in Fig.|S) and an upper limit for 
the extension along the line of sight, L„a^-, we can get an es- 
timate of the lower limits of «/ and «//. The shocked inhomo- 
geneities cannot extend beyond the Vela SNR, so the length of 
the chord intercepted by the Vela shell along the line of sight 
represents the maximum extension L„,a v of the clouds in this 
direction. If we approximate the Vela shell as a sphere with ra- 



dius rveia «i 5.06 X 10''' cm, in FilD L, 



5.8 X 10''' cm and 
5.1 X 10'^ cm and 



nj > 0.32 cm~^, while, in RegNE, L„,a 
nil > 0.14 cm"^. 

On the other hand, we may assume that FilD has an average 
extension along the line of sight comparable to its size in the 
plane of the image (LfHD ~ 2x 10'^ cm) and, in this way, we 
obtain a best guess value for «/ ~ 1.5 cm"-^. In the same way, 
setting LRegNE ~ 3 x 10'^ cm, we obtain a best guess value 
for nil ~ 0.55 cm"^. Notice that these values indicate pressure 
equilibrium between cores and coronae. 

An independent estimate of a lower limit of «/ can be ob- 
tained from the evidence that the cloud cores have reached 
the collisional ionization equilibrium. Modelling the cooler 
compo nent of the spectra with the PSHOCK non-equilibrium 
model llBorkowski et al.l l2001 ), we obtained best-fit values of 
T]siEi ~ 2 X lO'^ s/cm~^ with a lower limit (for subregion 4) of 
2 X 10" s/cm"^ so: 



tnei 



nh'ela 



ni(t)dt = niitveia - tFiio) > 2 X 10" s/cm^ (1) 



where {tveia - tpiiD) is the elapsed time since the shock reached 
FilD, and 77/ is the average density in this period. Assuming the 
Sedov model describes the evolution of the shock in the inter- 
cloud medium, we have: 



tVela - tpilD 



r''"° dr 



tVela 



1 



rvda 



5/2 



(2) 



If FilD and the center of the Vela shell lie on a same plane, 
normal to the line of sight, rfHD » 4.1 x 10'^ cm , therefore 
tveia - tpiiD ~ 4500 yr and, from Eq. O, «/ > 1.4 cm""*. 
Conversely, by putting «/ = 5 cm""' we can also obtain a lower 
limit for tveia - tpuD, in this way we infer that FilD was shocked 
at least 1300 yr ago. 

A similar useful estimate of a lower limit for nu was not 
possible, because in the coronae we do not have a univocal sig- 
nature of CIE and the hotter component can also be described 
with a NEI model, as we show in Sect. 14.61 Our results can be 
summarized as follows: 



Cores: 



«/ « 1.5 (0.32 H- 5) 
H/ ^ 1.4 (cm^^) 



(cm-3) 



Coronae: nu 0.55 (0.14 h- 5) 



(cm^^) 



4.3. 3D map of the plasma 

Merging the results of the spectral analysis with the related 
images, we derived the three-dimensional structure of the ob- 
served ISM clouds with resolution down to less than 0.1 pc. 
These maps are based on the assumption that the observed X- 
ray emission at any position in the field of view can be de- 
scribed with two-component thermal models having the same 
temperatures values, but different emission measure weights. 
The two temperatures were computed as mean values over all 
the spectral subregions presented in Sect. 13.31 also assuming 
the same Ne/Fe ratio anywhere. This approach allows us to de- 
rive maps of the extensions along the line of sight of the cores 
(L/) and of the coronae (L//) 
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In each pixel, the MOS count rates Rs, in the 0.3 - 0.5 keV 
band, and /?/„ between 0.5 keV and 1 keV are: 



4-7T 



Rh = Pt 



- n]LiQ. 



'(Ti) 



4n 



+ PZ 



An 



An 



(3) 



(4) 



where P^I^^{T) is the function of radiative losses per emis- 
sion measure folded with the XMM-Newton MOS response (in 
the corresponding energy band) and Q. is the solid angle sub- 
tended by a pixel. The P'^^^ (T) values were calculated with the 
simulation of a MOS observation (in the soft and in the hard 
band) of an isothermal plasma with temperature T///, abun- 
dances Fe/FeQ and Ne/Ne^ and with {rPjLiQ.IAn) = 1. In this 
way we found: Pf^^ij,) = 20.08 cm^ s"'; Plf^Hfi) = 14.21 
cm^ s-i; Pfo^(f//) = 53.53 cm^ s"'; andPf^^(T//) = 191.02 
cm^ s"'. 

Solving the system of Eqs. (|3} and (0} for jjLjjj we ob- 



tam: 



(An 



n^Li = I — 



RhPf°HTn)-RsPrHTn) 



=- (5) 



2 /47!- 



RsPr' 



■RhPf'^'iTi) 



pMos (Tu)Pr' (T,)- Pr' (T,)Pr' (Tu) 



■(6) 



Notice that Eqs. (|3} and © hold only if the count-rate in each 
pixel does not depend on the count rate in the adjacent pixels, 
so the angular extension of the pixels has to be larger than that 
of the PSF. For this reason we realized two MOS count-rate im- 
ages (in the 0.3 - 0.5 keV and 0.5 - 1 keV bands respectively) 
by rebinning with a pixel-size of 16". We then applied Eqs. (|5ji 
and (|6j at each pixel of these images and thereby obtained a 
high resolution map of n^L for the two components. From this 
map it is possible to know the emission measure of the plasma 
with a spatial resolution % 0.05 pc. It is important to remark 
that the results obtained from Eqs. (|5} and are consistent 
with the findings of the spectral analysis; for example, in sub- 
region 14 (where Tu is significantly different from Tu) the best 
fit value of n]jLu is 1.13+0.05 x lO'** cm"^, while from the rp-L 
map we obtain that, in this subregion, the mean value of n^,L// 
is - 1.09x10"* cm-'*. 

Assuming our best values for the density of the cloud cores 
and of the cloud coronae, «/ = 1.5 cm"^ and «// - 0.55 cm"^ 
(Sect. 14.2'1 . respectively, we derived the extension along the line 
of sight of the two components for each pixel of the rP'L map; 
Figs. 1^ and 1101 show the 3D structure of the cores (in white) 
and of the coronae (in transparent blue), presented from two 
different angle-shots. 

In the map the large core of FilD is surrounded by a very 
thin corona, the two peaks correspond to the two bright parts of 
the FilD (A and B in Fig.[2and subregions 2 and 4 of Fig.|3J, 
and the optical filament of Fig|3llies between them. In RegNE 
the core is nearly absent and the corona is very thick. It is possi- 
ble that the core of RegNE lies outside of the EPIC FOV, but we 
can also suppose that this cloud is just behind the front shock 



and that its core has not yet been reached by the shock (see also 
Sect. 14. 6> . Compared to FilD, the southwestern structure has a 
smaller core and a thicker corona. 

We can also derive an estimate of the mass of FilD. The 
dimensions of its core are ~ 2.4 x 10'^ cm and ~ 1.1 x 10'^ cm 
(on the plane of the observation) and 6.8 x 10'**/«^ cm (along 
the line of sight), while its corona has about the same extension 
on the plane of the observation; but it is an order of magnitude 
less thick, so the mass of FilD, for «/ - 1.5 cm"^, is ~ 10^' 
g (the mass of the corona is ~ 7 x 10^^ g). The dimensions of 
RegNE are - 1.4xl0'^cmand ~ 3.3 x 10'^ cm, on the plane of 
the observation, and ~ 1 x lO'^/n^^ cm, along the line of sight. 

Note that with other guesses of the «/ value, assuming pres- 
sure equilibrium, the extensions along the line of sight of cores 
and coronae scale with the same factor (L oc 1/n^); hence, the 
relative size of the inhomogeneities does not change. 

It is possible to show that the 3D map is consistent with the 
hypothesis of pressure equilibrium between cores and coronae. 
In fact the time-scale for pressure balance is the sonic time- 
scale. Therefore we expect pressure equilibrium into structures 
with dimensions comparable to the length A covered at the 
sonic speed Cj = {jkbT/fiy^^ in 4500 yr, which is the time 
elapsed from the shock-impact. A/ x 1.8 x 10'^ cm for the 
cooler component, and A// » 2.8 x 10'^ cm for the hotter 
component; and these values are indeed comparable to the to- 
tal (core-Hcorona) extensions of the observed inhomogeneities, 
thus supporting the hypothesis of pressure equilibrium. 

4.4. Shock-cloud interaction 

In this section we discuss the dynamics and evolution of the 
shock-cloud interaction in the observed region in light of our 
new results. 

The results of our analysis allow us to conclude that, since 
we neither see enhancement of the X-ray emission "behind" 
the optical filament nor inhomogeneities in the plasma temper- 
ature and density, the reflected shock model (see Sect.[2i is not 
compatible with our data. This result is further supported by 
the trend of the points in the E vs. R scatter-plot, as explained 
in Sect.lO 



4.4.1 . Cloud crushing time 

An important time-scale for evoluti on of the shock-clo ud inter- 
action is the cloud crushing time tec llKlein etalJl990l) . defined 
for a spherical cloud with radius p as 



PX 



1/2 



(7) 



where, for FilD, p x 1.2 x 10'^ cm and;if is the ratio of the 
density of the cloud to that of the inter-cloud medium. Since 
Vjc/x^^^ is the velocity of the shock into the cloud, this is the 
time-scale for the transmitted shock to cross the cloud. 

To obtain an estimate of x, we found an upper limit for the 
inter-cloud density This value was obtained introducing a 
third hotter MEKAL component in the analysis of the spectra 
of the subregions 12 and 16, where there are no distinct signa- 
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tures of the presence of relevant inhomogeneities, and consid- 
ering the upper limit of its emission measure. In this way we 
found n,e $ 0.06 cm"^. Since the density of FilD is ~ 1.5 
(Sect.^,X ~ 30. 

We can estimate v,c using the Sedov model, in this case: 



v(ic) 



^' Vela -ill 
jtVela 



(8) 



where rpuD > 4.1 x 10''' cm (Sect.O and rpm < 5.0 x lO''' 
cm. This upper limit is obtained from Eq. Q with tveia - (fud ~ 
1300 yr and then from Eq. ^. Therefore v,, = 5.7 - 7.8 x 10^ 
cm/s. In this way we get 

2700 < tac < 3700 yr. 

We have already pointed out (Sect. 14. 2> that FilD was 
shocked less than 4500 yr ago. Moreover, assuming nj =1.5 
cm""*, we obtain that the time elapsed from the shock impact is 
about 4200 yr; hence, we conclude that we are observing the 
shock-cloud interaction at ^ Itcc- 



4.4.2. Thermal conduction and radiative cooling 

The evolution of the shocked ISM clouds, apart from dragging, 
is governed by competition of radiative cooling and thermal 
conduction with the surrounding medium, so it is important to 
estimate the efficiency of both these processes to describe the 
scenario of the shock-cloud interaction. 

The radiative cooling time-scale of an optically thin plasma 
is defined as: 



Tr = 



(r - l)n.n//P(r) 



(9) 



where y is the the specific heat ratio. The va lues of P(T) in the 
0.01 - 50 keV band (calculated accord ing to lMewe etalJl985[ 
iMewe etalJl 19861 iLiedahl et al]ll995h for the temperature of 
the cores and of the coronae and for the observed abundances, 
are pCT;) « 1.0 x 10 "^^ and P(r,/) ^ 3.4 x lO"^^ gj-g ^-i cm\ 
Considering the 'S pitzeJ ( 1 19621) thermal conduction, the 
thermal conduction time-scale is: 



Ink 



1 



2(r-i)^.r5/2(i//2. 



(10) 



where /,■ indicates the temperature scale-height in the three di- 
mensions and ks = 5.6 x 10"^ erg s"' cm"' K"^^^ Since in the 
observed inhomogeneities we detected only two temperatures, 
we can reasonably approximate the temperature scale-heights 
with the extensions of the clouds. The dimensions of FilD and 
of RegNE are reported in Sect. 14.31 the southwestern structure 
extends beyond the EPIC FOV, so it is not possible to estimate 
its dimensions. 

We estimated and Tc for the cores and the coronae of 
FilD and of RegNE for all the permitted values of the density 
(Sect. 14. 2> . assuming pressure equilibrium between cores and 
coronae. The coronae are heated through thermal conduction 
from the hotter and less dense inter-cloud medium; on the other 
hand, they lose their thermal energy through radiative losses 
and thermal conduction towards the cooler cores. Notice that. 



since the inter-cloud medium is hotter and less dense than the 
coronae, Eq. ilO\ predicts that the heating of the coronae is 
more efficient than their cooling through thermal conduction. 

In FilD, for all the permitted values of density, we found 
Tr w 4 X 10^ yr and t^. ^ 2 x 10^ yr: the radiative cooling is 
negligible, so the FilD corona gains thermal energy by means 
of thermal conduction from the inter-cloud medium. This heat- 
ing produces the progressive evaporation of the corona. It is 
possib le to estim ate the rate of evaporation of the corona using 
the McKee & Cowie.(J977i) relation for a spherical cloud with 
radius ppc (in pc) : 



dm 
dt 



2.75xl0^rfp;,,(30/lnA) g/s 



(11) 



where m is the mass of the cloud and In A is the Coulomb log- 
arithm (for r > 4.2 X IQ^K, In A = 29.7 + ln(«:^'^^(7W10*')). 
Assuming T,,. = 6 x 10*" K, n,^ = 0.05 cm"^ (Sect lTTlT i and 
Ppc - 0.5 X (/f -I- ly + l,)l?>, we obtain an evaporation time-scale 
Tev - m/{dm/dt) of ~ 100 yr. This means that the FilD corona 
has almost totally evaporated in the few thousand years elapsed 
from the shock impact, in agreement with Fig. |5] and Fig. 1101 
where this corona is thinner than the others. 

The FilD core is heated through thermal conduction from 
the surrounding hotter corona. On the other hand, it loses ther- 
mal energy ffirough radiative losses and, if we assume that 
the optical filament Ues inside FilD, through thermal conduc- 
tion toward the optical filament. The radiative and conductive 
cooling time scale are invariably of the same order (t, /t"c = 



0.8 - 
and 



1.5) and are ~ 5 x 10 yr (~ 



10 yr, for «/ - 1.4 cm 



2 X 10 yr, for «/ = 5 cm ). The heating time-scale 



is shorter and corresponds to the thermal conduction cooling 
time-scale of ffie FilD corona (~ 2000 yr). It is possible to esti- 
mate whether this core is gaining or losing energy: the heating 
rate is // = (3nj[kTijVjj)/Tc, where Vu is the volume of the 
FilD corona and t^ is ffie conductive time-scale of the corona; 
the cooling rate is C (3njkTiVj)/Tr, where Vj is the volume 
of the core and is its radiative cooling time-scale. The ratio 
C/H is always lower than 0.22, so we conclude that ffiere is a 
net heating of the FilD core. 

Since the characteristic time-scale of energy exchanges be- 
tween the core and corona are comparable with ffie sonic time- 
scale, the hypothesis of pressure equilibrium is self-consistent. 

As for RegNE, where the core is nearly absent, for all the 
permitted values of density, ~ 500tc. Hence the radiative 
cooling is also negligible in this cloud. Using Eq. il 1> we esti- 
mated that ffie time necessary for ffie evaporation of ffiis corona 
to be completed is ~ 500 yr. 

4.5. The physical scenario in FilD 

As we pointed out in the previous sections, the thermal conduc- 
tion between the cores and the coronae and especially between 
the coronae and the inter-cloud medium is a very efficient pro- 
cess. On the other hand, the estimated value of tec (Sect. imi 
and the morphology of ffie optical filament within the FilD core 
(and into the southwestern structure) indicate ffiat a transmitted 
shock is propagating through a cloud with an increasing den- 
sity profile CBocchino et al...2000.) . The emerging scenario is 
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that the observed X-ray emission comes from the cloud plasma 
heated "internally" by the transmitted shock and "externally" 
by the thermal conduction from the hotter inter-cloud medium. 
In the inner and denser part of the cloud, the shock is radia- 
tive and the plasma cools off through radiative optical emis- 
sion, while in the cores and in the coronae there is a net heating 
of the plasma which evaporates in the inter-cloud medium in a 
few 10^ yr. Fig.lTTI is a schematic representation of the struc- 
ture of the FilD cloud, which summarizes all the ISM phases 
we identified on the basis of optical and X-rays data; we also 
report the values of temperature and particle density of each 
phase. 

4.6. NEI effects 

In all the extracted spectra, the fittings of the cooler compo- 
nent with the PSHOCK non-equilibrium of ionization model 
jBorkowski et al.li200ll) revealed that the cores have reached 
the collisional ionization equilibrium, because the best-fit val- 
ues of the tnei parameter are a few lO'-' s/cnT^ (and their lower 
limits are always larger than a few lO" s/cm~^), and the best- 
fit temperatures and emission measures are the same as those 
obtained with the MEKAL model. 

The statistics of the hotter component are too low to in- 
fer significant information on NEI effects in the FilD and in 
the southwestern coronae, while in RegNE (subregions 14 and 
15) the spectra are described well by the non-equilibrium of 
ionization model with a temperature of ~ 7 x 10'' K. The re- 
sults obtained in this region with the CIE and NEI conditions 
are reported in Table|2] Assuming NEI conditions for the hotter 
component, the temperature of RegNE is a bit higher, and com- 
parable to the inter-cloud temperature (Bocchino et al. 1999). 
This is at odds with the patchy morphology of the emission 
of this region, which seems to suggest that RegNE cannot be 
associated with the inter-cloud medium. 



Table 2. Best-fit parameters of the spectral analysis of subre- 
gions 14 and 15 with the MEKAL CIE model and with the 
PSHOCK NEI model. 



MEKAL 




subregion 14 


subregion 15 


Temperature (10*' 


K) 


2.58+^™ 


9 7^+0.13 


n-L{lO" 




1 1 3+"^ 


1 1 2+"-7 






64+" ''' 


54+" " 


xHd. 0. f.) 




146.4(125) 


158.7(134) 


PSHOCK 




subregion 14 


subregion 15 


Temperature (10^ 


K) 


' -2 


7!^ 


n^L(10" cm-') 




2 3+' 


2.5«, 






49+" 07 


"■-'"-0.10 


Tnei (s/cm^) 




6!* X 10'" 


8+5 X 10'" 


X'(d. 0. f.) 




131.2(124) 


145.6(133) 



A more detailed study of NEI conditions will require more 
data; however, we stress that our conclusions are not affected 
by NEI. In fact, since the temperature of the hotter component 
is higher than in CIE, the thermal conduction to the cores is 
more efficient, so the time-scale for the heating of the cores is 
lower and the evaporation rate of the coronae is higher than in 
CIE. 

However, from preliminary analysis of the archive EPIC 
observation 0137550901, pointing the Vela SNR at coordinates 
a = 8''48'"54' and 6 = -45°39'37", we found that in that re- 
gion of the Vela shell the plasma can be described with two 
thermal components at ~ 1 x 10^ K and 3 x 10* K, respectively, 
in equilibrium of ionization. This result indicates that the cloud 
coronae can be in equilibrium of ionization. 

5. Summary and conclusions 

Analyzing the EPIC observation of the Vela FilD region, we 
derived a detailed description of the morphology of the emit- 
ting plasma. In particular, the data allowed us to link the optical 
and X-ray emission of the FilD cloud. The optical filament lies 
between the two FilD regions with the highest X-ray surface 
brightness and this correlation, which seems to be present also 
in the southwestern structure, indicates that the emission orig- 
inates from ISM inhomogeneities with an inward increasing 
density profile. 

The study of the relationship between the observed count 
rate and the mean photon energy allowed us to separate the dif- 
ferent inhomogeneities present in the FOV. The spectral anal- 
ysis, performed on the spectra extracted from 16 subregions, 
showed that the plasma is described well by an optically-thin 
plasma in CIE with two thermal components. The two com- 
ponents are associated with two different phases of the clouds: 
the cloud core (which corresponds to the cooler component), 
surrounded by a hotter and less dense corona. The temperature 
of the cloud cores is uniform in the FOV (~ 1.14 + 0.02 x 10*" 
K) and also all the cloud coronae have the same temperature 
(3.0 + 0.1 X 10* K). The similarity of the temperatures implies 
almost uniform densities among the cores and among the coro- 
nae. The inhomogeneous X-ray emission is instead determined 
by large variations of the emission measure of both compo- 
nents. This indicates different volume distributions along the 
line of sight and allowed us to derive a 3D distribution of the 
observed plasma for the first time. 

The observed X-ray emission cannot be associated with 
the inter-cloud rn edium, and the reflected shock model 
llHesteret al.ll 19941) is not compatible with our data. The ISM 
inhomogeneities are heated by the transmitted shock which, 
travelling through the clouds, generates X-ray and optical emis- 
sion. 

We also addressed the role of thermal conduction and radia- 
tive losses in the post-shock evolution of the clouds. We found 
that thermal conduction is much more efficient than radiative 
losses, and it causes heating of the clouds. This heating deter- 
mines the evaporation of the outer layers of the clouds; after 
about one cloud crushing time the corona of FilD is almost 
wholly evaporated. 
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We also derived the mean abundances, between all subre- 
gions, of O (O/O0 * 1.0), Fe (F^IFe^ = 0.39 + 0.05), and Ne 
(Ne/NeQ = 1.7 + 0.2). The Fe abundance can hardly be recon- 
ciled with the post-shock dust grain destruction model. We do 
not have a straightforward explanation for Ne overabundance. 

These results allow us to deepen our knowledge of the in- 
teraction of an SNR shock with a small 1 pc) isolated inho- 
mogeneity. It is then interesting to verify if our conclusions can 
also be generalized for other regions of the Vela shell. Other 
authors have pointed out that physical conditions of the Vela 
SNR environments are very different, e. g. the southwestern 
part v s. the northeastern part. For instance, iNichols & SlavinI 
( I2OO4I) noted that the presence of large dense clouds can ex- 
plain their discovery of a band of strong high excitation C 1 near 
the western front of the Vela SNR, where the X-ray emission 
decreases and many optical filaments are present. Mor eover, 
the ^•^C O survey of the Vela SNR presented by Moriguch i et alJ 
ll200lh shows that the large scale spatial correlation between 
the X-ray emission and the molecular clouds seems to be quite 
complicated, while iLu & Aschenbachi tOOd have shown that 
the temperatures and EMs of the two thermal components nec- 
essary to describe the X-ray spectra present significant inhomo- 
geneities across the remnant. We then expect that the physical 
parameters we found in FilD may not be representative of the 
entire Vela SNR. In this framework, in a forthcoming paper, we 
aim to perform a systematic study of small regions, selected on 
the basis of large surveys, to shed light on the correlation be- 
tween X-ray emitting plasma and other ISM components (dust, 
CO , etc.). 

Future studies should also address other open problems, 
like the lack of NEl effects and Fe and Ne abundances. 
Moreover, to better understand the physics of the shock-cloud 
interaction, our further investigations will be focussed on the 
comparison between the observation results and those that will 
be obtained by performing detailed hydrodynamical simula- 
tions. 
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Figure 1. Adaptively smoothed count rate images (MOS equiv- 
alent) in the 0.3 - 0.5 keV band (left) and 0.5 - 1 keV band 
(right). These images are a weighted average of the pn, MOSl, 
and MOS2 images. Eight contour levels, equispaced between 
0% and 100% (included) of the peak value (~ 7.1 x 10"^ cnt 
s"' bin"' and ~ 10"** cnt s"' bin"' respectively), are shown; the 
bin is 4". 



Figure 2. Left: Adaptively smoothed count rate images (MOS 
equivalent) in the 0.3 - 2 keV band (bin size= 4"); Eight con- 
tour levels, equispaced between 0% and 100% (included) of the 
peak value (~ 1.7 x lO"'* cnt s"' bin"') are shown. The white 
square indicates the field of view of Fig. Right: Mean energy 
map (smoothed with cr - 20", bin size= 10"). In each pixel the 
local mean photon energy of the photons detected by the MOS 
cameras between 0.3 keV and 2 keV is reported. We superim- 
posed 6 contour levels equispaced between 0.5 keV and 0.85 
keV. The 16 subregions selected for spectral analysis are also 
shown. 



Figure 3. Color com posite optical image of FilD (from 
iBocchino e t al."2000): in the electronic version, the Ha emis- 
sion is reported in green, and the [Olll] emission in violet. The 
5.7' X 5.7' field of view is indicated by the white square in the 
left panel of Fig. |2l We superimposed (in red) X-ray contour 
levels derived from EPIC observation in the 0.3 - 05 ke V band. 



Figure 4. E vs. count rate scatter plot for the pixel of the left 
panel of Fig.|2l In the electronic version the colors of the points 
codify the three different regions indicated in the field of view. 
Errors in £ ~ 10 eV, bin size= 10". 



Figures, pn (upper) and MOS-summed (lower) spectra ex- 
tracted from subregion 7 with the corresponding best fit model 
and residuals. 



Figure 6. Best fit Ne abundances in units of the solar photo- 
spheric abundance (the errors are at 90% confidence level) ver- 
sus the mean MOS count rate (in bin of 10") for the 16 sub- 
regions of Fig. 12] The Fe abundances can be deduced from 
the condition (Ne/NeQ)/(Fe/FeQ) - 4.4. CrosssFilD sub- 
regions; diamondsssouthwestern subregions; squaressRegNE 
and North subregions. 
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Figure 7. Best fit values of Tj (cooler component) e Tn (hotter 
component) versus the mean MOS count rate (in bin of 10") 
for the 16 subregions of Fig.|2| In subregions 14 and 15 (whose 
spectra are well described by one component) only Tu is indi- 
cated. T errors are at 90% confidence level, while the maximum 
errors in the count rate are of 3 x 10"* s"'. Symbols as in Fig. 

m 

Figure 8. Best fit values of the emission measure per unit area 
(that is, the product of the square of the density n of the plasma 
times its geometric extension along the line of sight L) of the 
cooler component, left panel, and of the hotter component, 
right panel, vs. the mean MOS count rate (in bin of 10") for 
the 16 subregions of Fig.|2l For subregions 14 and 15 (whose 
spectra are well described only by the hotter component), an 
upper limit to the cooler component is reported. Errors as in0 
symbols as in Fig. |6l FilD = subregions 1-6; S-W = subre- 
gions 7 - 12; RegNE and North region = subregions 14 - 16; 
the subregion 13 is between FilD and RegNE. 



Figure 9. 3D map of the observed ISM clouds. The height of 
the features is proportional to the extension along the line of 
sight of the cores (L/, in white) and of the coronae (L//, semi- 
transparent blue). North and East are also indicated. In this im- 
age the structure of FilD is visible; the optical filament of Fig. 
lilies between the two peaks, corresponding to subregions 2 
and 4 of Fig.|2l Assuming «/ = 1.5 cm"-' and «// = 0.55 cm"-', 
the maximum value for Lj is ~ 4.4 x 10'^ cm and for Lu is 
~ 4.7 X 10"^ cm (for other values of density L scale as 1 /«^). 

Figure 10. Same as Fig.|9lfrom a different angle-shot. The large 
corona of RegNE is visible in close-up. 

Figure 11. Schematic representation of the structure of the FilD 
cloud. We report the values of the particle density n and of the 
temperature T of the plasma; the value s of the optic al filament 
temperature and density are taken from lBocchino et aL (t200Qi) . 
The fat arrows indicate the evaporation of the plasma. 

Appendix A: Influence of the background on the 
spectral best-fit parameters 

As explained in Sect. 12. 21 we used the EPIC public observation 
0106260201 to extract a pn and a MOS background spectrum 
for the spectral analysis. 

In Table IaI we report the best fit parameters obtained from 
the spectral analysis of subregions 4, 10, and 16, with the sub- 
traction from source spectra of the background of the EPIC 
observation 0106260201 and of the EPIC observations of the 
Galactic Plane Survey'. The table shows that the choice of the 
photon background does not significatively influence the values 
of the best-fit parameters. 



Table A.l. Best fit parameters, obtained from the spectral anal- 
ysis of subregions 4, 10, and 16 of Fig. |2] with subtraction 
from source spectra of the background of the EPIC observa- 
tion 0106260201 (left column) and of the background of the 
10 EPIC observations of the Galactic Plane Survey (GPS, right 
column). 



subregion 4 


0106260201 


GPS 


T, (10^ K) 


1.14+"-!!^ 

-0.04 


1.15+"-!!* 

-^-0.05 


T,i (10* K) 


3.4 + 0.5 


3.5+"-^ 


n]Li (10"* cm"^) 






niLii (10"* cm"^) 

JJ 11 \ " / 


0.28+"-!;:! 

-0.07 


0.28+"-?'' 

-0.1 


Ne/Neo 




1 3+' 


X' (d- o. f.) 


102.1 (85) 


93.0 (85) 


subregion 10 


0106260201 


GPS 


Ti (10'' K) 


1 05 +!;-!;^ 


1 01 +"-"„' 


T,i (10* K) 


-^-0.3 


3 1+"--' 

-0.1 


)i^L; (10"* cm"^) 


•^-0.7 


5.1+'-^ 


n\,L,i (10"* cm"^) 


0.52+»-!il 

-0.04 


0.54+"-"^ 

-0.06 




1 9+0.4 
^■^-0.3 


1 4+0.4 


(d. o. f.) 


181.0(157) 


176.9 (157) 


subregion 16 


0106260201 


GPS 


Ti (10* K) 


1.3+0.2 




Tu (10* K) 


2.9 + 0.2 




n\Li (10"* cm"^) 


16+"-* 

"•^ -0.03 


19+"-'^ 


4L„ (10"* cm-5) 


0.14 + 0.01 


14+"-"' 

"■^^-0.01 


NelNe^ 


^■"-0.2 


9 9+0.5 
^■^-0.4 


X' (d. 0. f.) 


189.1 (169) 


185.5 (169) 



' 10 EPIC observations of the Galactic Plane Survey, pointing 
galactic longitude 319°, and galactic latitude between 0° and 3°. 



